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Abstract 

Absorption and emission spectra of plasmas with multicharged-ions 
contain transition arrays with a huge number of coalescent electric- 
dipole (El) hnes, which are well suited for treatment by the unresolved 
transition array and derivative methods. But, some transition arrays 
show detailed features whose description requires diagonalization of 
the Hamiltonian matrix. We developed a hybrid opacity code, called 
SCORCG, which combines statistical approaches with fine-structure 
calculations consistently. Data required for the computation of de- 
tailed transition arrays (atomic configurations and atomic radial in- 
tegrals) are calculated by the super-configuration code SCO {Super- 
Configuration Opacity)^ which provides an accurate description of the 
plasma screening effects on the wave-functions. Level energies as well 
as position and strength of spectral lines are computed by an adapted 
RCG routine of R. D. Cowan. The resulting code provides opacities for 
hot plasmas and can handle mid-Z elements. The code is also a power- 
ful tool for the interpretation of recent laser and .^-pinch experimental 
spectra, as well as for validation of statistical methods. 



1 Introduction 

Computation of plasma opacities is necessary in many fields of plasma 
physics: experiments are diagnosed with emission or transmission spectra 
and mean opacities are required to evaluate radiative transfer. Both ap- 
plications require calculation of spectra with a high resolving power, i.e., 
of over a few hundreds, which involves computation of atomic structure 
and transition arrays. However, plasma spectroscopy rarely requires one to 
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get over a resolving power higher than a few thousands, thus fine structure 
modehng provides sufficiently accurate spectra. 

Nowadays, plasma opacity is calculated by codes which use DLA (De- 
tailed Line Accounting) like HULLAC [1] and FAC [2] or, for faster compu- 
tation, statistical methods based on UTA (Unresolved Transition Array) or 
STA (Super Transition Array) approaches. 

DLA approach gives precise enough transition rates and energies for 
spectroscopy and mean opacities calculation, but requires substantial re- 
sources as soon as configurations get more complex (degeneracy gc ^ fOOO). 
However, the DLA approach alone enables one to account for the effects of 
the Boltzmann factor e'^^^-' on probability of levels 7 J within a configura- 
tion, even if statistical approaches can still be improved [Gilleron et al., this 
issue]. The DLA approach is also useful when transition arrays have small 
total physical broadening, due to various mechanisms, e.g., Doppler, Stark, 
collisions. 

On the other hand, statistical approaches enable much faster opacity 
calculations, because order to 2 moments of transition arrays can be ob- 
tained without diagonalization of the Hamiltonian matrix. That is the basis 
of UTA (non-relativistic) [3j and SOSA (relativistic Spin-Orbit Split Array) 
[1] approaches. UTA and SOSA still require detailed configuration account- 
ing (DCA). In case of a large number of configurations, the STA approach 
[5] enables an even faster calculation of plasma opacities, replacing config- 
urations by superconfigurations. SCO [6] is able to handle both STA and 
DCA for opacity calculation. However, statistical approaches are limited by 
the spectral accuracy they can reach, and, thus, are not able to resolve all 
structures experimental spectra can show. Further, statistical approaches 
neglect the effects of temperature on probability distribution over levels of a 
single configuration. This probability distribution can have a strong impact 
on transition array shape if the plasma temperature is of the same order of 
magnitude as the configuration energy span, or receptive zone. 

However, it is worth investigating the relevance of a DLA computation 
of very complex transition arrays. As shown in the examples in table [U 
complex transition arrays with millions of lines can occur in mid-Z plasmas 
at moderate temperatures, e.g., copper at 20 eV, for doubly-excited configu- 
rations, when the Sd subshell is almost half-filled. In cases when a / subshell 
is partially opened, even singly-excited configurations show a very complex 
structure. Since statistical and detailed approaches are both useful for an 
accurate description of spectra over a broad range of temperature, density 
and atomic number, a model mixing both approaches could substantially 
improve the quality of the opacity calculation, without requiring too many 
computational resources. 
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Ion 


Transition array 


Number of lines 


Fe XIV at 150 eV 
Cu IX at 40 eV 
Gd VI at 40 eV 


2p^3s^3p3d - 2p^3s^3p3d:^ 
2p^3p^3d^ - 2p^3p^3(f 
Afbd - 4p5d6d 


177,684 
564,293 
1,139,911 



Table 1: Number of spectral lines of transition arrays starting from singly- 
or doubly-excited configurations. 

2 The hybrid approach 

The hybrid approach we use in our model should be able to compute accurate 
spectra within a reasonable time, mixing configurations and super configura- 
tions at local thermodynamic equilibrium (LTE). Bound-bound opacity may 
include detailed transition arrays, as well as UTA, SOSA and STA. The cal- 
culation should be as consistent as possible from a physical point of view. 
Effects of the plasma environment on atomic structure should be considered 
even in DLA computations. We also intend to provide an automatic tool, 
able to select the (super)configurations to be studied as well as to decide 
how to handle each transition array. 

2.1 Implementation 

We choose to use two existing codes, SCO and R. D. Cowan's RCG routine 
[7]. SCO is a superconfiguration opacity code based on a physical model that 
gives a reliable description of plasma physics over a wide range of tempera- 
ture, density and atomic number. A special feature of SCO is the automatic 
selection of configurations and super configurations, starting from an average 
atom model and using fluctuation theory [8l |^. Atomic structure is com- 
puted using a self-consistent approach within central-field approximation, 
and density effects like free-electron screening are taken into account by the 
code. Orbital relaxation can be included for selected transitions. There are 
relativistic corrections based on Pauli approximation, which enables correct 
handling of medium to high Z elements. 

On the other hand. Cowan's code can compute detailed transition arrays, 
given an ion and its electronic configuration. Atomic structure and atomic 
spectra are calculated by two different programs, respectively RCN/RCN2 
and RCG. Output of the former is automatically formatted as input (atomic 
structure integrals) for the latter. The separation between the calculations 
of atomic structure and of transition arrays enables the input of atomic 
structure integrals computed by other processes, e.g., SCO. We choose to 
adapt RCG as a subroutine in SCO only used for DLA computation. 
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2.2 Description of SCORCG 



The resulting code is able to handle consistently both detailed and sta- 
tistical transition arrays. Input is a Fortran 77 list of variables in which 
thermodynamic quantities, a list of subshells needed, the number of super- 
configurations and the ionization range, among other options, can be set. 
Most parameters have a default setting suitable for standard calculations, 
thus input can be rather small. An average atom plasma is computed at 
first. It enables one to have an overlook at the atomic structure and other 
global properties of the plasma to be studied. For the computation itself, 
this first run selects the configurations that are to be included in the opac- 
ity calculation. The code can be stopped after this first run and the list of 
selected configurations can be kept as is or modified. 

In SCORCG, a superconfiguration E can match a single configuration. 
Such single superconfigurations are treated with the DCA methods. Other 
super configurations are specifically called non-DCA superconfigurations. For 
each superconfiguration S in the list generated as input, SCO processes 
self-consistent calculations of atomic structure, with relativistic corrections, 
free-electron screening inside ions. Effects of neighbouring ions are limited 
by the Wigner-Seitz sphere. Orbital relaxation can be included if wished 
[lOj . Next, SCO computes all transition arrays starting from H. 

Only transition arrays starting from single configurations can be detailed. 
SCO handles all cases of UTA, SOSA and STA, and sends to the modified 
RCG routine arrays that are to be treated in detail, according to criteria 
driven by the physical broadening as well as computational issues. The 
spectral opacity of the elemental plasma can be written as the sum of the 
opacities of the species in it, weighted by their respective probabilities: 



where H is any superconfiguration in the plasma, its probability, and 

Ksihiy) its spectral opacity was H the only species in the plasma. 

The modified RCG routine does not only compute line strengths and 
energies with radial integrals given by SCO for the detailed configuration 
C, it also calculates other features of the detailed transition arrays, such as 
strength- weighted moments of the distribution of lines. Moments of orders 
up to 100 can be calculated in statistical weight approximation (SWAP) like 
in Eqs. ([8]) and ([9]) and by using the Boltzmann factor. The detailed parti- 
tion function Uc of the starting configuration is also calculated according to 
Eq. ([2]) by our modified RCG routine. Each transition array starting from 
C is then multiplied by Uc- 




(1) 
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2.3 Management of different types of species in plasmas 

The general expression for the partition function of the configuration C is 

Uc= Y.i2J + l)e-^^''\ (2) 
-rJeC 

where 7J is a fine-structure level of C, characterized by its total angular 
momentum J and its energy E^j; (3 is the inverse of temperature kT. This 
partition function is used in SCORCG when at least one transition array 
starting from the configuration is calculated in detail. 

If all transition arrays starting from C are statistically treated, then C 
will not be detailed. Statistical treatment of configurations and transition 
arrays assumes the probability of each level 7J does not depend on tempera- 
ture. Therefore, the partition function of C only depends on the population 
Ws of the subshell s characterized by its angular momentum ig-. 

Uc=Y. (2-^ + ^y^'''' = 5ce-^^^ (3) 

with 

7Jec s ^ ^ 

After the computation of the opacities of all superconfigurations, the to- 
tal bound-free and bound-bound opacities are divided by the total partition 
function lA of the plasma: 

C detailed C statistical H non— DCA 

Hence, the opacity of each superconfiguration H in Eq. ([1]) is multiplied 
by the probability of the latter, was H detailed, DCA or non-DCA: 

= (6) 

2.4 Current limitations of the model 

There are two types of limitations in SCORCG: the first are computational 
while the second are due to the model. Little computational resources are 
required by a STA or UTA code, and one can run such opacity programs on 
a desktop computer. DLA greatly increases the complexity and slows down 
computation if the detailed configurations are complex. 

The computational limitations are propagated to the generation of the 
opacity from our modified RCG routine. Hence the DLA approach is limited 
to transition arrays with less than 800,000 lines, a limit that can be set lower 
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upon input file. For one transition array, the size of blocks with the same J 
in Hamiltonian matrix is limited to 4,000. However, this limitation is much 
less restrictive than the restriction of lines in a transition array. Moreover, 
a maximum of 8 subshells can be opened in both configurations. 

The plasma is modeled using average-atom and (super) configurations. 
Fine-structure information is only used for calculation of lines or partition 
functions; atomic structure and modeling of plasma density effects are statis- 
tically treated. Relativistic effects are accounted in the Pauli approximation, 
which provides accurate results for mid-Z elements. Full configuration inter- 
action (CI) is not implemented in SCORCG. Only configuration interactions 
between relativistic subconfigurations that belong to the same configuration 
are handled in the calculation (exactly in DLA part and by an approximate 
method in the statistical part). 

DLA calculations are much more time-consuming than UTA or STA. For 
most transition arrays, a great deal of computational resource is consumed 
in the convolution of transitions with physical line profiles. The first solution 
we implemented is to regroup all lines of a transition array into energy bins 
in which spectral line intensities are added. In SCORCG, we assume that 
all lines of the same transition array have the same broadening as the statis- 
tical transition array for which SCO routines compute Stark and collisional 
formulas [11]. Even if "binning" reduces the amount of resources required 
for a single transition array with many lines, numerous lines remain to be 
computed within a whole opacity calculation. For further improvement of 
DLA efficiency, one can use one average broadening per bin for all transition 
arrays. 

Abdallah et al. [12] showed this approach gives moderately accurate 
spectra, in that Rosseland mean opacity changed by about 1 or 2 % from 
a calculation with a detailed line by line convolution. However, the method 
they developed for two broadening bins per energy bin - one for narrow 
lines and one for wide lines - gives much more accurate results, with almost 
no observable difference for the opacity spectrum and yielding a Rosseland 
error of about 0.01 %. The latter method requires one to sample all spectral 
lines twice: first for computation of average broadenings for each energy bin, 
and second for filling bins with spectral lines that fit the criteria of energy 
and broadening. 

3 Experiment interpretations 

We have analyzed experimental spectra produced in laser and Z-pinch fa- 
cilities, focusing on spectra for which there is no satisfactory interpretation 
done when statistical codes were used. Three spectra for which our hy- 
brid model strongly improves the agreement with experimental data are 
presented. The first spectrum was recorded in an experiment at the HE- 
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LEN laser facility in the United Kingdom [13]. The second was obtained in 
an experiment at Sandia Z-pinch facility in the USA[T^. The third, most 
recent, is a transmission spectrum recorded at LULI 2000 laser facility of 
Ecole Polytechnique in France jl5]. 

3.1 Aluminum on HELEN laser facility 

The aluminum experiment on HELEN laser facility was thought to be at 
a temperature of 40 eV and a density of 0.01 g/cc [T3]. The experimental 
transmission shows spectrally-resolved features that were not well repro- 
duced by any statistical approach. Interpretation of this spectrum using 
SCORCG gives a much better agreement, with a temperature of approx- 
imately 40 eV if density effects are not accounted, and of 36 eV if these 
effects are added. Although not understood, the differences we observe 
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Figure 1: Aluminum spectrum from the experiment of Davidson ^13j com- 
pared with SCO and SCORCG calculations. 

3.2 Iron on Z-pinch 

The iron spectrum produced at Sandia National Laboratory (SNL) [T^ has 
a high resolution and shows very detailed features. This transmission spec- 
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trum was intended to provide a benchmark spectrum which could test opac- 
ity simulations at a temperatures above 100 eV and electron densities higher 
than 10^^ cm~^. This spectrum shows features no statistical approach is able 
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Figure 2: Iron spectrum of Bailey [14j and the best agreement found with 
SCO and SCORCG. 

to reproduce. Only codes using a DLA approach can obtain good agreement 
with the experimental spectrum, and the results already obtained by other 
codes , e.g., PrismSPECT [TB] make interpretation challenging. The best 
agreement for the SCORCG calculation with experiment is at T = 150 eV 
and p = 58 mg/cm^. Some features do not seem to be reproduced by our 
mixed detailed and statistical approach, but ion by ion investigations show 
the discrepancies do not come from ionic distribution. 



3.3 Copper on LULI 2000 laser facility 

For the LULI experiments we mainly focus on shot 31 of 2008 campaign for 
opacities of mid-Z elements |15] . Shot 31 had a copper sample heated to 
16 eV and with a density of 5 mg/cc. The larger areal mass (40 ;Ug/cm'^) 
makes this transmission spectrum more difficult to model, as bound-bound 
structures are prevalent. 

Even if SCO gives good agreement, the ratio between relativistic sub- 
structures 2pi/2 ~ 3^3/2 and 2^3/2 — 3^5/2 at 12.8 and 13.1 Angstroms is not 
reproduced by the fully statistical model. But the SCORCG transmission, 
like other DLA codes [Blenski et al., this issue], gives good agreement with 
the experiment. It is an effect of relativistic configuration interaction be- 
tween the two major SOSAs of 2p—3d transition. However, a new relativistic 
CI model ^17J shows a much better agreement of SCO with experiment and 
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Figure 3: Transmission spectrum of copper obtained in shot 31 [TC] and 
theoretical spectra computed by SCO and SCORCG. 

DLA codes. 

4 Array shapes and moments 

With its ability to compute transition arrays in both statistical and DLA 
approaches, SCORCG allows us to perform detailed investigations of the first 
to fourth order moments of the transition arrays and their consequences for 
the spectra. 

4.1 Effects of order 3 and 4 moments on modeling of transi- 
tion arrays 

Line distribution in transition arrays is modeled by a Gaussian distribution 
in any UTA approach. However, if spin-orbit splitting effects are significant, 
the array is split into several subarrays. Even when spin-orbit splitting 
is very low, the distribution of lines in the transition arrays is often not 
symmetric and neither the central peak nor the wings of the transition array 
are well-fit by a Gaussian, as shown in Fig. HI We define the order n moment 
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of a transition array with this formula: 



7J-7'J' 

/^n = . (7) 

/ , ^-yJ-^f'J' 

^j-yj' 

S^j-yj' is the area of the transition array and /ui is the average 

7J-7'J' 

energy. For order 2 and higher, we define centered moments as: 

{Eyj-yj' - fJ-iT ■ S^j-yj' 
l^^' = ^ , • (8) 

■yj-yj' 

(c) 

H2 is the variance of the transition array. For higher orders, we use the 
reduced centered moments of the transition array, i.e., 

Order 3 moment is the skewness of the array, which represents the asym- 
metry of the hne distribution. Order 4 moment is the kurtosis of the array, 
which represents the sharpness/flattening of the distribution near its cen- 
ter. Pain et al. [18] demonstrated how third and fourth orders moments 
can significantly improve the modeling of transition arrays in a statistical 
approach. Figure H] shows that inclusion of the asymmetry and kurtosis im- 
proves the modeling of transition arrays with low spin-orbit splitting. The 
normal inverse Gaussian, NIG, agrees much better with the shape of the 
wings of the transition array, and this improved profile can have an impact 
on Rosseland mean opacity calculated by statistical approaches. 

Figure [S] shows that most transition arrays, even those that contribute 
most to the opacity, are asymmetrical and have a broad range of values 
for the kurtosis. Transition arrays seem to converge toward non-Gaussian 
shapes when n grows. Those starting from subshell Is tend to concentrate 
near (as, 04) = (+0.8,2.5). 

Accounting for the order 3 and 4 moments can strongly improve the 
statistical modeling of some transition arrays with a very peaked shape. 
These moments can be obtained by a DLA calculation, knowing that for 
transition arrays with a large number of lines, a large part of computation 
time is spent in convolving all lines by Voigt profiles. 



10 
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Figure 4: Transition array 3p^3dP — Sp^Sd^Ap from Ge XII with broadening 
of 0.1 eV and temperature of 50 eV. Comparison DLA/UTA-Gaussian/NIG. 
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Figure 5: Distribution of order 3 and 4 moments used in the interpreta- 
tion of Davidson's aluminum experiment. Areas enclosed by circles, squares 
and triangles are proportional to the product of the total transition array 
intensity and the abundance of the initial configuration. 
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4.2 Effects of temperature on average energy and variance 



Temperature, or more exactly the inverse of temperature (3 



kT 



has a 



significant effect on shape of transition arrays. Gilleron et al. in Ref. [this 
issue] show significant impact of temperature on the total absorption of a 
transition array, even at temperatures higher than the statistical width of 
the transition array. Investigations into effects of temperature on the average 
energy and dispersion of the transition array show that the average energy 
does not seem to depend much on temperature, because all lines are in a 
limited energy range. But the spectral resolution required for an accurate 
spectral calculation is about 1000, whereas temperature impact on average 
energy of the transition array is about 10 % for some transition arrays. 
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Figure 6: Variances of selected transition arrays in Fe XII. 



Figure [6] shows that dependence is quite strong for variance, which can 
be partly explained by the fact that at low temperature, only transitions 
starting from ground level of the initial configuration are significant. This 
lowers the dispersion of the transition array. The moments of the transition 
array - 3p^3d, of type £^~^^ — £'^£', are not much affected by temperature. 
However, the moments of transition arrays starting from configurations with 
two open subshells are significantly affected by low temperature. This can 
be explained by a wider span in energy for the starting configuration. 
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4.3 Effects of the maximum number of lines per transition 
array on mean opacities 

As a harmonic mean, the Rosseland mean opacity is sensitive to the pres- 
ence and the treatment of spectral hnes when physical broadening does 
not merge them into smoother structures. With our new hybrid code, we 
started to investigate the effect on the Rosseland mean opacity of changing 
the DLA treatment. We focused on one interesting case: copper at 20 eV 
and 10~^ g/cc. At this temperature, the Rosseland mean opacity is ex- 
tremely sensitive to the amount of DLA treatment of the transition arrays. 
Between a calculation with as much DLA as possible and a fully statistical 
calculation, the Rosseland mean differs by 20 %, as shows Fig. [71 




Maximum number of lines per ti'ansition aiTay 



Figure 7: Rosseland mean opacity plotted against the maximum number of 
lines per transition array 

Figure [3 shows that the Rosseland mean opacity only converges to «;/j(oo) = 
2.135 cm^/g, the limit of mean opacity with every transition array detailed, 
if transition arrays with more than 100,000 lines are treated in detail. Al- 
though kr(oo) seems to be reached below 10,000 lines per transition array, 
some structures of the bound-bound spectrum, especially in the 3p — 3d 
region at around 60 eV are not converged at this number of lines per array. 

Spectral envelopes in Fig. [5]explain the peaks between 25,000 and 65,000 
lines per transition array. In comparison, a similar test on an aluminum 
plasma shows that convergence of Rosseland mean opacity occurs at a few 
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maximum 800,000 lines 

■ - ■ maximum 50,000 lines 

■ ■ ■ ■ maximum 10,000 lines 
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Figure 8: Envelopes of spectra of plasma opacities with different maximum 
number of lines (upper bounds merge). 

hundreds of lines per transition array. 

5 Conclusion 

SCORCG, the hybrid code we are developing, is a versatile tool for inves- 
tigations on plasma opacity computation. Interpretation of experimental 
spectra has given results that are comparable to those obtained with other 
DLA codes. Beyond this first application, SCORCG enables us to study 
more deeply statistical properties of bound-bound transitions in hot plas- 
mas by the collection of substantial data, which allows us to automate the 
execution and produce spectra in a reasonable time. Therefore, SCORCG 
will be quite useful in a future investigation into the validity of statistical 
approximations for opacities in hot plasmas, as it enables one to see their 
effects on the entire spectrum. However, this code is still under develop- 
ment, and further improvements involve spectral line handling, criteria for 
discrimination between transition arrays that require DLA calculation or 
accounting of order 3 and 4 moments. 



14 



References 

[1] A. Bar-Shalom, M. Klapisch, and J. Oreg, J. Quant. Spectroc. Radiat. 
Transfer 71, 169 (2001). 

[2] M. F. Gu, Can. J. Phys. 86, 675 (2008). 

[3] C. Bauche-Arnoult, J. Bauche, and M. Klapisch, Phys. Rev. A 20, 2424 
(1979). 

[4] C. Bauche-Arnoult, J. Bauche, and M. Klapisch, Phys. Rev. A 31, 2248 
(1985). 

[5] A. Bar-Shalom, J. Orcg, W. H. Goldstein, D. Shvarts, and A. Zigler, 
Phys. Rev. A 40, 3183 (1989). 

[6] T. Blenski, A. Grimaldi, and F. Perrot, J. Quant. Spectrosc. Radiat. 
Transfer 65, 91 (2000). 

[7] R. D. Cowan, The Theory of Atomic Structure and Spectra, The Re- 
gents of the University of California, 1981. 

[8] F. Perrot, Physica A 150, 357 (1988). 

[9] B. G. Wilson, J. Quant. Spectrosc. Radiat. Transfer 49, 241 (1993). 

[10] F. Perrot, T. Blenski, and A. Grimaldi, J. Quant. Spectroc. Radiat. 
Transfer 58, 845 (1997). 

[11] M. S. Dimitrijevic and N. Konjevic, Astron. Astrophys. 172, 345 (1987). 

[12] J. Abdallah. Jr., D. Kilcrease, N. Magee, S. Mazevet, P. Hakel, and M. 
Sherrill, High Energy Density Phys. 3, 309 (2007). 

[13] S. J. Davidson, J. M. Foster, C. C. Smith, K. A. Warburton, and S. J. 
Rose, App. Phys. Lett. 52, 847 (1987). 

[14] J. E. Bailey, G. A. Rochau, C. A. Iglesias, J. Abdallah, J. J. MacFarlane, 
I. Golovkin, P. Wang, R. C. Mancini, P. W. Lake, T. C. Moore, M. 
Bump, O. Garcia, and S. Mazevet, Phys. Rev. Lett. 99, 265002 (2007). 

[15] G. Loisel, P. Arnault, S. Bastiani-Ceccotti, T. Blenski, T. Caillaud, J. 
Fariaut, W. Flsner, F. Gilleron, J.-C. Pain, M. Poiricr, C. Reverdin, 
V. Silvcrt, F. Thais, S. Turck- Chieze, and B. Villette, High Energy 
Density Phys. 5, 173 (2009). 

[16] J. J. Mcfarlane, I. Golovkin, P. Woodruff, and P. Wang, APS Meeting 
Abstracts J, 1181 (2003). 



15 



[17] F. Gilleron, J. Bauche, and C. Bauche-Arnoult, J. Phys. B: At. Mol. 
Opt. Phys. 40 (2007). 

[18] J.-Ch. Pain, F. Gilleron, J. Bauche, and C. Bauche-Arnoult, High En- 
ergy Density Phys. 5, 294 (2009). 



16 



